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Abstract 

There are claims in the literature that in neutrino oscillations and oscillations 
of neutral kaons and U-mesons the oscillation phase differs from the standard one 
by a factor of two. We reconsider the arguments leading to this extra factor and 
investigate, in particular, the non-relativistic regime. We actually find that the 
very same arguments lead to an ambiguous phase and that the extra factor of 
two is a special case. We demonstrate that the unitarity triangle (UT) fit in the 
Standard Model with three families is a suitable means to discriminate between the 
standard oscillation phase and the phase with an extra factor of two. If Kl — K$ 
and B<iH — B^l mass differences are extracted from the K° — K° and f?9 — data, 
respectively, with the extra factor of two in the oscillation phases, then the UT fit 
becomes significantly worse in comparison with the standard fit and the extra factor 
of two is disfavoured by the existing data at the level of more than 3 a. 
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1 Introduction 



Compelling evidence in favour of neutrino oscillations obtained in recent years in the 
Super-Kamiokande [1,2], SNO [3], KamLAND [4], K2K [5] and other neutrino experiments 
(see e.g. [6] and references therein) is a major breakthrough in the search for physics 
beyond the Standard Model. All existing neutrino oscillation data with the exception 
of the LSND data [7] 1 are well described if we assume three-neutrino mixing. Defining 
^ m |fc = m i — m h where the rrij are the neutrino masses, the best fit values 

Am 2 21 = 7.9 x 10~ 5 eV 2 and | Am 2 2 | = 2.4 x 10~ 3 eV 2 , (1) 

were found for the solar [4] and atmospheric neutrino neutrino mass-squared differences [9] , 
respectively. 

These values of the neutrino mass-squared differences were obtained from neutrino 
oscillation data under the assumption that the neutrino transition and survival probabil- 
ities have the standard form (see e.g. the reviews in Ref. [10]). Neutrino oscillations are 
due to the interference of the amplitudes of the propagation of neutrinos with different 
masses and the standard phase differences are given by the expression 

&Pj* = -jf-- ( 2 ) 

Here E is the neutrino energy and L is the distance between neutrino production and 
neutrino interaction points. The theory of neutrino oscillations has a long history start- 
ing with the paper of Gribov and Pontecorvo [11] (for other early papers see [12,13], for 
historical overviews see [14]). There is also a rich literature on more elaborate derivations 
of neutrino transition and survival probabilities based on quantum mechanics and quan- 
tum field theory (for a choice of these papers see [15-25], more citations are found in the 
reviews [26-29]), which all result in the standard oscillation phases of Eq. (J2J). 

There exist, however, claims [30] that the phase differences in neutrino transition 
probabilities differ from the standard ones by a factor of two and are equal to 

Am%L 



A<p jk = (3) 

Other authors [31] claim that there is an ambiguity in the oscillation phase. Theoretical 
discussions about the factor of two or other factors in oscillation phases continue during 
many years — see e.g. [23,25,29,32,33] where these additional factors have been refuted on 
theoretical grounds. Taking into account the fundamental importance of the problem we 
believe that it is worthwhile to think about possibilities to confront the different oscillation 
phases to experimental data. 

The same non-quantum-theoretical arguments which lead to an additional factor of 
two in neutrino oscillation phases can be applied to the oscillation phases in M° ^ M° 
oscillations of neutral bosons M° = K°, B®, etc., as was demonstrated in Ref. [33]. A 
more complicated additional factor has been obtained in Ref. [34], but was subsequently 

1 The result of the LSND experiment is planned to be checked by the MiniBooNE experiment [8] which 
is currently taking data. 
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refuted in Ref. [35]. Since in M° M° oscillation experiments the mesons are often 
non-relativistic, the relevant oscillation phase is 

where p is the momentum of the neutral meson. In the ultra-relativistic limit, Eq. (j^J) 
coincides with Eq. (J2J). In the following we use the subscript QT for the standard phase (@J), 
whereas phases different from the standard phase are marked by a bar — see Eq. (JHJ). 

In recent years a remarkable progress in the measurement of \V c b\, \V u b\ and other 
elements of the CKM matrix was reached (see e.g. [36]). Another great achievement 
was the measurement of the CP parameter sin 2(3 with an accuracy of about 5% in the 
BaBar [37] and Belle [38] experiments at asymmetric B-factories. This allowed to perform 
a new check of the Standard Model based on the test of the unitarity of the CKM mixing 
matrix, the so-called unitarity triangle test of the SM. It was shown [39-43] that the SM 
with three families of quarks is in an good agreement with existing data, which include the 
data on the measurements of the effects of CP violation. In the unitarity triangle (UT) 
test the experimental values of the Ki — K$ mass difference Amx and the B^u — B^l 
mass difference Am^ are used. The values of Amx and Ams d were obtained from an 
analysis of the experimental data based on the standard transition probabilities with the 
standard oscillation phase (jlj). 

In this paper we will present the result of the UT test under the assumption that 
oscillation phases in K° ^=> K° and B^ ^ oscillations differ from the standard ones 
by the above factor of two. We will show that such an assumption is disfavoured by the 
existing data at the level of more than 3 a. 

The plan of the paper is as follows. In Section|2]we will discuss in some detail how this 
notorious factor of two in the oscillation phase appears. Considerations how to confront 
the factor of two with experiment are found in Section |3J Section 0] contains our UT fit 
with and without the factor of two. Our conclusions are presented in Section |S| The 
technical details of the UT fit are deferred to an appendix. 



2 The notorious factor of two 
2.1 Notation 

For simplicity we consider oscillations between only two states. Thus we have two different 
masses rrij (j = 1,2). We adopt the convention mi < m 2 . For each mass eigenstate the 
relevant phase is 

<Pj = E i t ~ Pj L i ( 5 ) 

where Ej = Jpj + m| and pj are energy and momentum, respectively. Though there are 
some arguments that in particle oscillations mass eigenstates with the same energies are 
coherent [20,21,25,33], we want to be general and assume neither equal energies nor equal 
momenta. 
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It is useful to define quantities Ap and Am via 



Pi,2 =PT ^Ap, mi^ = m =f ^Am, (6) 

where p and m denote average momentum and mass, respectively. Defining Am 2 = 
m 2 , — mf and Am = 1712 — mi, we have the relation 

Am 2 = 2mAm. (7) 

In the following we will use the approximations 

p ^> |Ap| with Ap = aAm. (8) 

The dimensionless constant a is zero for pi = p 2 . In general it will be of order one or 
even larger. In the non-relativistic case one can have a ~ m/p. The first relation of 
Eq. (JHJ) excludes particles which are nearly at rest; such a situation is not contained in 
our discussion. Consequently, we do not allow for p « m or a » 1. However, we will 
take care that all our considerations hold also in the moderately non-relativistic limit. 
The second relation in Eq. (JHJ) states our coherence assumption: mass eigenstates with 
momenta which differ more than the mass difference can be coherent. Note that with 
Eq. (jHJ) we have 

p > Am. (9) 

In the following we will need 

! Am? uAt) I 

AE = E 2 -E 1 = - (mAm + pAp) = — — + ^-f with E = - {E 1 + E 2 ) . (10) 
E IE E 1 

2.2 "Derivation" of extra factors in oscillation phases 

Particle oscillation phases different from that of Eq. (J3J) have been found for instance 
in Refs. [30,34], and an ambiguity of a factor of two in the oscillation phase has been 
diagnosed in Ref. [31]. It was stressed first in Ref. [33] and then in Refs. [29,32] that in 
essence the discrepancy to the standard result (J1J) is due to the assumption that the two 
mass eigenstates are detected at the same space point but at different times 

tj = L/vj = LEj/pj. (11) 

For each mass eigenstate, the corresponding time tj is inserted into the phase (0). The 
motivation for this is that particles with different masses move with different velocities 
Vj. This picture mixes quantum-theoretical and classical considerations in an ad hoc 
fashion and leads to the conclusion that particle phases taken at different times, though 
at the same space point, produce the interference, which is in contradiction to the rules 
of quantum theory. 

Eq. (|TT|) gives the phase 

Pj Pj 
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and, therefore, the phase difference 



P2 Pi 

Then, using only Ap <C p, we obtain 

A^2A PQT - M±|W . ( 14 ) 

As seen from this equation, Aip differs from Ay?Q T not only by a factor of two, but also 
by an additional term which contains the arbitrary quantity 2 Ap. In the ultra-relativistic 
case, which always applies to neutrinos but also to M° £5 M° oscillations when their 
energy is high enough, the additional term is negligible and we have the ultra-relativistic 
phase 

(A^) UR ^2A^ QT . (15) 

For oscillations of non-relativistic neutral flavoured mesons, the additional term can not 
only be comparable with the first term but could even dominate in Eq. (fT^j) . Since Ap 
is arbitrary, we come to the conclusion that, for oscillations of non-relativistic particles, 
Eq. (jllj) leads to an arbitrary — and thus unphysical — oscillation phase. 

In order to illustrate the latter point, let us consider the two extreme cases of equal 
momenta and equal energies. In the first case with Ap = 0, Eq. (|14j) gives 

- — Am 2 L 2mAmL , . 

A<p = = . 16 

p p 

Clearly, we have again the notorious factor of two, in comparison with the quantum- 
theoretical result. On the other hand, equal energies correspond to Ap = —Am 2 / (2p) 
(see Eq. dTUJ)) and with Eq. (JUI) the result is 

—r — Am 2 L ( ' m 2 \ 

A^ = 1 + 7TT • ( 17 ) 



P V 2p 2 

This oscillation phase, which is similar to the one advocated in Ref. [34], agrees with 
Eq. (|16|) only in the ultra-relativistic limit. 



2.3 The quantum- theoretical oscillation phase 

Although it has been stressed many times (see e.g. Ref. [27]) that the quantum-theoretical 
oscillation phase does not suffer from any ambiguity, it is instructive to repeat the deriva- 
tion of this fact here, in order to compare with the derivation of Eq. (|14|) . Quantum 
theory requires the two phases (0) to be taken at the same space-time point. Therefore, 
we have 

A</3qt = AET — ApL, (18) 

In principle, one should be able to determine an upper limit on Ap from the widths of the wave packets 
of the particles participating in the neutrino, K°, B^, etc. production and detection processes [18,20,24]. 
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where T characterizes the time when the interference takes place. Then, with T = LE/p 
we obtain the quantum-theoretical result 



/ Am 2 pAp\ EL Am 2 L mAmL 

A ^ = {^ + E w)T- ApL = — = — ' (19) 

where the arbitrary quantity Ap has dropped out. 3 For M° t> M° oscillations, the 
phase (|T9~j) can also be written in the familiar form Am r, where r is the eigentime of the 
particle for covering a distance L. 

We want to emphasize that a more complete understanding of the oscillation phase 
needs a full quantum-mechanical or quantum field-theoretical approach. All such treat- 
ments (see for instance the reviews [26,28,29] and references therein) consistently give the 
result of Eq. (fTTJj) . In approaches not guided by quantum mechanics or quantum field the- 
ory the conversion of time into a distance is always the subtle point [33,35]. In all present 
experiments, oscillations are treated as phenomena in space. If eigentimes are used for 
the evaluation of data, then distances are converted into times (see e.g. [37,38,44]). 



3 Confronting non-quantum-theoretical phases with 
experiment 

Since we have seen that the derivation of phase (pPfj) does not conform to the rules of 
quantum theory whereas Eq. (jlj) does, then one could ask the question why consider the 
phase (fTl^ at all. From our point of view, the reason for this is twofold: 

• On the one hand, there is the subtlety that the time difference At = \t 2 — ti| (see 
Eq. (HU)), which is the culprit of the discrepancy with the quantum-theoretical 
result, is immeasurably small. 

• On the other hand, as we will show, the phases (fTB^) and (fT7|) can actually be tested 
experimentally. 



The time difference can be expressed as 

L 



At 



2pE 



Am 2 — (m\ + wij) — — 



(20) 



To get a feeling for the size of At, we take the K°K° system with Am^ ~ 3.48 x 10~ 12 MeV 
and use for example L — 1 m, p — 1 GeV and Ap = 0. Then we find At ~ 5 x 10~ 24 sec, 
which is indeed far beyond measurability. 

As for an experimental test of the phase (fTTj) we consider two different measurements 
of the Kl — Kg mass difference. Since this phase has an additional dependence on the 
momentum, it is useful to compare two measurements which have different average kaon 



3 It is reasonable to assume that T is Ljv\ or L/v2 or some average of these two expressions. What 
one takes precisely as T is irrelevant, because all these possibilities differ only in terms suppressed by 
Am and Ap. Since AE is already small in that sense (see Eq. HlOjl) and the first order in Am and Ap is 
sufficient, we take the velocity p/E. 
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momenta. The CPLEAR experiment has measured [45] Arn^ = (5295±20±3) x 10 6 Hs -1 . 
In that experiment kaons are produced in the reaction pp — » K + tt~K° and the charged- 
conjugate reaction, with pp annihilation at rest. Thus the kaons are non-relativistic. In 
the KTeV experiment the kaons are in the ultra-relativistic regime; this experiment has 
obtained [46] Attlk = (5261 ± 15) x 10 6 hs^ 1 . According to Eq. (|17jl the mass differences 
extracted in these experiments should be different and related by 

(Am^) CPLEAR _ ^ m 2 R0 y i + m K° (21) 
(Am K ) KTeV 2p 2 KO ~ 2p 2 K0m J 

where px° is the (average) neutral- kaon momentum in the CPLEAR experiment. 4 One 
can show that the maximal energy of the neutral kaon in the CPLEAR reaction is given 
by 

4?Hp- m 2 — 2 m n m K 

^°max = : , (22) 

4m p 

where m p , m w and mx are proton, pion and kaon mass, respectively. For our purpose the 
distinction between the mass values of the charged and neutral kaon masses is irrelevant. 
With the numbers above for the mass differences obtained by the CPLEAR and KTeV 
experiments, we use the law of propagation of errors to compute the value 1.006±0.005 for 
the ratio on the left-hand side of Eq. (|2ip. We insert the values of the particle masses into 
Eq. (J22J and calculate Pxo max ; then we arrive at 1.22 for the right-hand side of Eq. (}2~T]) . 
which is about 40 standard deviations larger than the ratio of Ki — K$ mass differences. 
Consequently, we conclude that the phase (JT7J) is in contradiction to the results of the 
CPLEAR and KTeV experiments. 

The phase (fTSj) which contains the notorious factor of two needs a different approach; 
in the next section we will use the fit to the unitarity triangle constructed from the CKM 
matrix to show that this factor of two is experimentally strongly disfavoured. For the idea 
to compare the Am 2 result of the solar neutrino experiments with that of the KamLAND 
experiment see Ref. [47]. 



4 The unitarity triangle fit 

4.1 Description of the unitarity triangle analysis 

Following the traditional way, the unitarity triangle (UT) is given by the three points 
A = (p,fj), B = (1,0), C = (0,0) in the plane of the parameters p and fj, which are 
defined by 

?=p(i-y). a = .i(i-y). (23) 

where A, p, rj are the Wolfenstein parameters of the CKM matrix. Pedagogical introduc- 
tions to the UT can be found e.g. in Refs. [39-41]. Our numerical analysis is based on 
the input data as given in Tab. 1 of Ref. [43], and we use the following constraints to 
determine the point A = (p, fj) : 

4 If Eq. (|17fl were correct, there should also be a dependence of the extracted mass difference on Pk°- 
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The measured value of Ek = (2.280 ± 0.013) x 10 3 . The theoretical prediction for 
this quantity, which is a measure for CP violation in K° — K° mixing, is given by 5 



(24) 



where Amx is the Kl — K$ mass difference and Bk, C,D,E are numbers which 
have to be calculated and/or depend on measured quantities such as A, m t , m c , \V c b\ 
(see e.g. Ref. [40] for precise definitions). 

The experimental determination of IK^/V^I. This ratio is connected to p, fj by 



(25) 



• The measurement of the B^u — BdL mass difference 

Am Bd = 0.502 ± 0.006 ftps" 1 . (26) 

The theoretical prediction for the square root of Am^ as a function of p, fj is given 
by 

v/A^; = F \V cb \ Xy/fj 2 + (1 - P) 2 , (27) 

where F is a constant depending on mt and other quantities subject to theoretical 
uncertainties (see e.g. Ref. [40]). 

• In addition we use direct information on the angles of the unitarity triangle a, j3, 7. 
The angle /3 has been measured at BaBar and Belle, and we use the value sin 2(3 = 
0.726±0.037. For 7 we use the value (59.1±16.7)° (see Ref. [43] for details), whereas 
for a we use the likelihood function extracted from Fig. 10 of Ref. [43] to take into 
account the two allowed regions for a around 107° and 176°. 

We do not use the constraint from Am^ which usually is included in UT fits. The reason 
is that at present only a lower bound exists on Am^ , and therefore no further constraint 
is obtained for the oscillation phase with the extra factor of two. However, we remark 
that once an upper bound on ArriB s will have been established in the future, this will 
provide additional information on the oscillation phase. 

The fit is performed by constructing a x 2- f unc tion x 2 {PiV) from these observables, 
including experimental as well as theoretical errors. Technical details on our analysis are 
given in Appendix EI 



4.2 Results of the UT analysis 

The result of the standard UT fit is shown in the upper panel of Fig. ^ It is in 
good agreement with the results of various groups performing this analysis, compare 
e.g. Refs. [36,43,48]. We show the allowed regions in the plane of p and fj at 95% CL for 

5 For the sake of brevity we drop the phase factor exp(i7r/4) in ejc, since it plays no role in the following. 
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Figure 1: Unitarity triangle fit with Arrix and Ams d obtained from the standard oscil- 
lation phase (upper panel) and the oscillation phase with the extra factor of two (lower 
panel). The shaded regions correspond to the 95% CL regions (2 dof) obtained from 
the constraints given by Ek, \V u b/V c b\, Am Bd and sin2/3. In addition, constraints from the 
measurement of the angles a and 7 are used in the fit (not shown in the figure). The 
ellipses correspond to the 95% CL regions from all data combined, and the stars mark 
the best fit points. 

the individual constraints from e K , |T4b/K&|> Am Bd , sin 2(3, as well as the combined anal- 
ysis including in addition the information on a and 7. The 95% CL regions are obtained 
within the Gaussian approximation for 2 degrees of freedom (dof), i.e. they are given 
by contours of A% 2 = 5.99. For the best fit point of the combined analysis we obtain 
p = 0.237,f/ = 0.325 with the 95% CL allowed region shown as the ellipse in Fig. [H 
Assuming that the x 2 -minimum has as a ^-distribution with (6 — 2) dof our value of 
Xmin = 1-4 implies the excellent goodness of fit of 84%. 

Let us now discuss how an extra factor of two in the oscillation phase will affect the 
UT fit. If such a factor is present the mass differences inferred from particle-antiparticle 
oscillation experiments will be two times smaller. Therefore, whenever in the UT analysis 
a mass difference inferred from oscillations enters one has to use 

Am = r Am (28) 
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0.5 1 1.5 2 2.5 

r = Am / Am 

Figure 2: x 2 oi the unitarity triangle fit as a function of the parameter r defined in 
Eq. ([28)1 . For fixed r the y 2 is minimized with respect to p and fj. 

with r = 1/2, where Am is the value obtained with the standard oscillation phase, i.e. 
this is the value which is given by the Particle Data Group [48]. In the lower panel of 
Fig. ^ we show the result of the UT fit by using the extra factor of two in the oscillation 
phase. This factor affects two observables relevant for the UT fit. 

1. In the prediction for ek shown in Eq. (J2"4^ the experimental value for Aitlk is used. 
Since this value is obtained from K° t> K° oscillations, Aitlk has to be replaced 
by ArriK if there is an extra factor of two in the oscillation phase. This moves the 
hyperbola in the (p, fj) plane from Ek to the right, as visible in Fig. ^ 

2. The experimental value for ArriB d given in Eq. (}2T)j) has to be replaced by ArriB d , 
which is a factor of two smaller. Therefore, from Eq. (|2T|) it is clear that the radius 
of the circle in the (p, fj) plane from Am^ d is reduced by a factor \/2, as can be seen 
also in Fig. ^ 

The other constraints from IK^/V^I, sin 2/?, a and 7 are obtained from particle decays 
without involving any oscillation effect, and therefore they do not depend on the oscillation 
phase. One observes from Fig. ^ that the agreement of the individual constraints gets 
significantly worse using the extra factor of two. In particular, at 95% CL there is only a 
very marginal overlap of the intersection of the allowed regions from IK^/V^I and sin 2(3 
with the one from e k . The best fit point in the lower panel of Fig. [T]has Xmin = 13.8, 
which implies a goodness of fit of only 0.8%, assuming a ^-distribution for 4 dof. 

In Fig. El we show the x 2 minimized with respect to p and f] as a function of the 
parameter r given in Eq. (}2"%|) . Hence, r = 1 corresponds to the standard oscillation 
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3.9 
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17.6 
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11.9 


3.8 
2.0 
3.2 


m c = (1.2 ±0.2) GeV 


1.4 


11.9 


3.2 


constraints on a, 7 not used 


0.13 


9.6 


3.1 



Table 1: The Xmin f° r the standard oscillation phase (r = 1) and for the oscillation phase 
with the extra factor of two (r = 1/2) for variations of the input data (see text for details). 
The column "number of a" gives the number of standard deviations with which r = 1/2 
is disfavoured with respect to r = 1. 



phase, and r = 1/2 corresponds to the extra factor of two. From this figure one observes 
the remarkable feature that the best fit point occurs exactly at r = 1. In other words, even 
if the extra factor in the oscillation phase is treated as a free parameter to be determined 
by the fit, the data prefer the standard oscillation phase. For the value r = 1/2 we obtain 
a Ax 2 = 12.4 with respect to the best fit point, which corresponds to an exclusion at 3.5cr 
for 1 dof. We conclude that the extra factor of two in the oscillation phase is strongly 
disfavoured by the UT fit. 

4.3 Robustness of the UT analysis 

In this subsection we investigate the robustness of our conclusion with respect to variations 
of the input data for the UT fit. To this aim we show in Tab. Qthe results of our analysis 
by changing some of the numbers entering the UT fit. The line "standard analysis" in the 
table corresponds to the analysis described in the previous two subsections. In particular, 
exactly the input data given in Tab. 1 of Ref. [43] are used. 

First we have investigated how our analysis depends on the value for \V u b\. We show the 
results of the fit by using only the value from exclusive (|T^j,[( exc i)) or inclusive (|T4b|(inci)) 
decays, where the numbers are taken from Ref. [43]. Note that in our standard analy- 
sis both values are taken into account, as described in Appendix [X] We observe from 
the numbers given in Tab. Q]that for the relatively small value for \V u b\ from exclusive 
measurements the fit gets notably worse for r = 1/2. In contrast, for the relatively large 
value from inclusive measurements the fit gets worse for the standard oscillation phase 
(Xmin = 3-9), whereas for r = 1/2 the fit improves with respect to the standard analysis 
(Xmin = 7.8). The reason is that for large values of \ V u b\ the radius of the circle in the (p, fj) 
plane from \V u b/V c b\ becomes larger, which worsens the fit for r = 1, whereas for r = 1/2 
the agreement of the individual allowed regions becomes better. Note however, that even 
for iVublfmci) the goodness of fit for r = 1/2 is only 1%, and r = 1/2 is disfavoured with 
respect to r = 1 by 2cr. We have also performed the analysis by using the (inclusive and 
exclusive) averaged value IK^pdg) obtained by the PDG [48]. The fit using the extra 
factor of two is slightly improved with respect to our standard analysis, however r = 1/2 
can still be excluded at 3.2cx. 

Furthermore we have investigated how our result depends on the input value for the 
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charm quark mass m c . The value m c = (1.2 ± 0.2) GeV is adopted by the CKM- fitter 
group [36], in contrast to the value m c = (1.3±0.1) GeV from the UTfit Collaboration [43] 
used in our standard analysis. The mild improvement of the fit for r = 1/2 comes mainly 
from the larger error on m c , which leads to a slightly larger allowed region from Ek- 

In the last line of Tab. Qwe have removed the constraints for the angles a and 7 from 
the fit, i.e. we use only ek, \ V u b/V c b\, AniB d , sin 2/3. We observe that the direct constraints 
of a and 7 contribute 4.2 units of x 2 to the Xmin for r = 1/2. However, also without the 
constraints for a and 7 the extra factor of two in the oscillation phase is excluded by more 
than 3a. 

Finally let us comment on the the very small value of Xmm = 0-13 (for 2 dof), which 
we obtain without the constraints on a and 7 for the standard oscillation phase. In 
fact, the x 2 -minimum value of 1.4 in the standard analysis comes mainly from a. To 
include the information on this angle we are using the likelihood function from Fig. 10 of 
Ref. [43] (see Appendix [21) , which has two maxima around 107° and 176°. The maximum 
at 176° is slightly preferred, whereas the UT fit requires the other maximum. The very 
small x 2 -minimum value obtained without using the likelihood for a shows that the fit 
is dominated by rather large theoretical errors. Therefore, x 2 is significantly lower as 
expected just from statistics. The fact that even with these assumptions on theoretical 
errors the x 2 is large for r = 1/2 implies that the exclusion of the extra factor of two in 
the oscillation phase is rather robust. 

5 Conclusions 

In this paper we have reconsidered claims that the standard oscillation phase (|19J) has to 
be corrected by extra factors. We have focused on possible tests of these extra factors by 
using experimental data. The usual starting point to derive these non-quantum-theoretical 
expressions for the oscillation phase is Eq. ([lip, which says that mass eigenstates with 
different masses need different times to reach the spatial point where the interference of the 
amplitudes for the different mass eigenstates takes place. In this way we have derived the 
phase Aip of Eq. (|14jh The aim of our theoretical discussion was to consider both neutrino 
oscillations and oscillations of neutral flavoured mesons. For M° — M° oscillations, it was 
important to include the non-relativistic limit in our phase considerations. 
We have obtained the following results: 

1. The non-quantum-theoretical phase Atp of Eq. (jUj) becomes ambiguous in the non- 
relativistic case, because it contains a small but arbitrary momentum difference 
Ap. We have stressed that in the correct quantum-theoretical treatment, where the 
amplitudes interfere at the same time, this arbitrary term does not show up. 

2. If we adjust Ap in Eq. (fPfj) such that the mass eigenstates have the same energy, 
then a momentum-dependent extra factor appears in Atp — see Eq. (|17|). We have 
shown that this extra momentum dependence is in disagreement with measurements 
of the Kl — Ks mass difference at different kaon energies. 

3. If Ap = 0, the notorious factor of two appears in Atp — see Eq. (II 6|) . We have 
demonstrated that using K L —K S and B dH — B dL mass differences extracted from the 
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data with the extra factor of two in the K° ^=> K° and B® ^ B® oscillation phases, 
respectively, the unitarity triangle fit in the Standard Model becomes significantly 
worse compared to the fit with the standard mass differences. The phase with the 
extra factor of two is excluded at more than three standard deviations with respect 
to the standard phase. 

Concerning this last point, as an additional check, we have treated the extra factor in 
the oscillation phase as a free parameter r (see Eq. ((2HJ)) and considered x 2 as a function 
of r. It is remarkable that the minimum of \ 2 occurs nearly precisely at r = 1, which 
corresponds to the standard oscillation phase. This result can be regarded as a successful 
test of quantum theory. It is likely that in the future, with accumulated data used in 
the unitarity triangle fit, the exclusion of the extra factor of two will become even more 
significant. 

Acknowledgements: S.M.B. acknowledges the support by the Italian Program "Rientro 
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of this work was performed. He is also grateful to A.Yu. Smirnov for a useful discussion. 
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A Details of our UT fit procedure 

The fit of the UT is performed by adopting the following x 2 -function: 

X \p,fi,B K ,\V ub \) = Y,( x T P -^ Cd )S^(xf p -xf cd ) + xl 

+ X 2 sy AB K ) + X 2 sy A\Vu b \) (29) 
The final x 2 is obtained by minimizing Eq. (|29p with respect to Bk and \V U b\: 

X 2 (p, fj) = Min [ X 2 (P, fj, B K , iKtl); Bk, iK&l] • (30) 

In Eq. (ffijj) the indices i,j run over (ex, \V u b/V c b\, Am^, (3, 7) and Sij is the covariance 
matrix of these observables containing the experimental as well as theoretical uncertain- 
ties. It also takes into account correlations between the various observables induced by 
the experimental errors of parameters such as m t , A and |V^|, which are common to more 
than one observable. The term \ 2 a contains the information on the angle a, and is defined 
as Xa — ~ 2 In [£ (a) /Max £(«)], where C(a) is the likelihood function for a read off from 
Fig. 10 of Ref. [43]. 

For the treatment of theoretical uncertainties we follow the common practice in UT 
fits to split the error into a Gaussian part and into a "flat" part, which cannot be assigned 
a probabilistic interpretation [36,42,43]. For the parameter Bk relevant for ek one has 
Bk = 0.86 ± 0.06 ± 0.14, where the first error is Gaussian and the second is "flat". To 
include both errors in our fit we construct a likelihood function C(Bk) by convoluting a 
Gaussian distribution with width 0.06 with a flat distribution which is non-zero in the 
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interval [—0.14, +0.14] and zero outside. Then this likelihood is converted into a x 2 by 
Xgy St (B K ) = —2 \n[C(B k ) /Ma.x£(B k )] which is added to the total x 2 according to Eq. (J2HJ). 
The resulting x 2 is minimised for fixed p and fj with respect to Bk- 

The value of \V u b\ can be obtained from exclusive and inclusive decays, where the 
exclusive measurement suffers from theoretical uncertainties characterized by a "flat" 
error (see e.g. Tab. 1 of Ref. [43]). In our standard analysis we include both values by 
constructing a likelihood function £ (\V ub \) = £ exc i(|K&|) x AnddKbl), where £ exc i(|K&|) is 
obtained similar as in the case of Bk by folding a Gaussian and a flat distribution, whereas 
Anci(|Kib|) is just a Gaussian distribution. Finally, the term Xsyst(l^l) i n Eq. flUj) is 
obtained by Xsystd^l) = ~ 21n[£(|V^&|)/Max£(|V! u &|)]. The dependence of our results on 
the treatment of \V u b\ is discussed in Sec. 14.31 
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